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Abstract 



(N 
> 

' We study possible back-reaction and quantum gravity effects in the evaporation of black holes 

, which could be produced at the LHC through a modification of the Hawking emission. The 

corrections are phenomenologically taken into account by employing a modified relation between 

\^ , the black hole mass and temperature. The usual assumption that black holes explode around 

■ 1 TeV is also released, and the evaporation process is extended to (possibly much) smaller final 

masses. We show that these effects could be observable for black holes produced with a relatively 

^ • large mass and should therefore be taken into account when simulating micro-black hole events 

I , for the experiments planned at the LHC. 

Qh, 
(D 

1 Introduction 

^ ■ One of the most exciting feature of models with large extra dimensions [H [2] is the fact that the 

^ \ fundamental scale of gravity Mq could be as low as the electro-weak scale [Mq ~ 1 TeV) and 

micro-black holes may therefore be produced in our accelerators O [U O [6l [7] . Once a black hole 
has formed (and after possible transients) the Hawking radiation II Ij is expected to set off. The 
most common description of this effect is based on the canonical □ Planckian distribution for the 
emitted particles and on the consequence that the life-time of micro-black holes is so short that the 
decay can be viewed as sudden [7j. This standard picture has already been implemented in several 
numerical codes |181 [T9l [20l [2T] which let the black hole decay down to an arbitrary mass Mf of order 
Mq via the Hawking law and then explode into a small number of decay products. By running these 
numerical codes, one mostly aims at extracting information that will allow us to identify micro- 
black hole events in the planned experiments based, for example, on some experimental features of 
the decay products [22] or on some particular signatures of the decay [23j . 
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^Black holes have been studied in both compact [5^ and infinite warped extra dimensions (see also Ref. [10]). 

^From a theoretical point of view, one should employ the more consistent microcanonical description of black hole 
evaporation [12] which was first applied in the context of large extra dimensions in Refs. [131 1141 116[ . It was then 
shown that actual life-times can vary greatly depending on the model details [161 117] . 
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We would like to stress here that the issue of the end of the black hole evaporation remains 
an open question (see, e.g., Refs. [Ml US |26]) because we do not yet have a reliable theory of 
quantum gravity and, most of all, no experimental data from the quantum gravity regime. The 
singular behavior of the Hawking temperature as the black hole mass decreases to zero can in fact 
be simply viewed as a sign of the lack of predictability of perturbative approaches. Therefore, 
although the very detection of these objects would already be evidence that Mq is not the Planck 
energy (about 10^^ GeV) and that extra dimensions may indeed exist, observing the late stages 
of black hole evaporation (when the black hole mass M ~ Mq) could provide us with the kind 
of data we need to finally build the theory of quantum gravity. On a purely experimental side, 
this is also an important issue, since it corresponds to determine whether some features usually 
associated with the black hole decay are sensitive to the assumptions made about the late stages 
of the evaporation. This is equivalent to asking whether deviations from the Hawking law, induced 
by an underlying and still unknown theory of quantum gravity, can actually be detected. In some 
sense, this is analogous to the problem of trans-Planckian modes in cosmology, where one can ask 
whether the still unknown physics at Planck scales, which is certainly important in the early stages 
of the Universe, can affect the CMB spectrum we observe now. 

As we mentioned above, the ignorance of the late stages of the black hole evaporation is bypassed 
in the numerical codes by letting the micro-black hole decay into a few standard model particles 
as an arbitrary lower mass is reached (the possibility of ending the evaporation by leaving a stable 
remnant has also been recently considered |27l I28j). The purpose of this report is precisely to 
study if modifications of this standard description of the evaporation which could occur as the 
black hole mass becomes smaller than some arbitrary scale (M < Mq) can produce experimentally 
detectable features, although we shall not consider a specific experimental setup and leave for 
future developments the task of including the detector's sensitivity and geometry. From another 
perspective, one can also view this preliminary analysis as an estimate of the systematic errors to 
be associated with the results already appeared on this subject that usually rely on the standard 
picture we described above. 

For this purpose, we have developed several modifications to the Monte Carlo code CHARYB- 
DIS [18j which we shall describe in Section [3] (see also Ref. [IS]). We shall then report some 
interesting and perhaps unexpected (preliminary) results from a few test runs in Section HI We 
shall use units with c = h = \ and the Boltzmann constant A;b = 1. 



2 Black hole evaporation 

No-hair theorems of General Relativity guarantee that a black hole is characterized by its mass, 
charges and angular momentum only. The one parameter characterizing an uncharged, non-rotating 
black hole is thus its mass M. On considering solutions to the Einstein equations (or applying 
Gauss' theorem) in 4 + d dimensions, one can then derive the following relation between the mass 
and the horizon radius, 

where F is the usual Gamma function, and the temperature associated with the horizon is given by 

d+l , , 
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Once formed, the black hole begins to evolve. In the standard picture the evaporation process can 
be divided into three characteristic stages 

1. Balding phase: the black hole radiates away the multipole moments it has inherited from 
the initial configuration, and settles down in a hairless state. A certain fraction of the initial 
mass will also be lost in gravitational radiation. 

2. Evaporation phase: it starts with a spin down phase in which the Hawking radiation [llj 
carries away the angular momentum, after which it proceeds with the emission of thermally 
distributed quanta until the black hole reaches the Planck mass (replaced by the fundamental 
scale Mq in the models we are considering here). The radiation spectrum contains all the 
Standard Model particles, which are emitted on our brane, as well as gravitons, which are 
also emitted into the extra dimensions. It is in fact expected that most of the initial energy is 
emitted during this phase into Standard Model particles [6| (although this conclusion is still 
being debated, see e.g., Ref. p9]). 

3. Planck phase: once the black hole has reached a mass close to the effective Planck scale Mq, 
it falls into the regime of quantum gravity and predictions become increasingly difficult. It is 
generally assumed that the black hole will either completely decay into some last few Standard 
Model particles or a stable remnant be left which carries away the remaining energy [27j . 

In our approach we will consider possible modifications to the second and third phases. On 
the one hand, we will allow a modified Hawking phase by employing a different relation between 
the horizon radius and the temperature. The necessity of such a modification, suggested by the 
fact that the Hawking temperature diverges as the mass of the black hole goes to zero, has been 
discussed in several papers, e.g., Refs. [301 [H]. We will follow a phenomenological approach to 
this issue and model a general class of modifications by making use of the form given in Eq. (j3.2p 
below. On the other hand, we will look at the possibility that the evaporation may not end at the 
fundamental scale Mq (~ 1 TeV) but proceeds further until a lower (arbitrary) mass Mf has been 
reached. 



3 Quantum Gravity and Monte Carlo code 

Several Monte Carlo codes which simulate the production and decay of micro-black holes are now 
available (see, e.g. Refs. \18 \ 119 1 [20]). We have found it convenient to implement our modifications 
into CHARYBDIS [18j. 

In order to let the black hole evaporate below the electro- weak scale Mewj at which the Standard 
Model particles acquire their mass, one needs to treat the mass of the emitted particles properly. 
In the original code, the phase space was taken to be that of massless particles since the mass of 
the black hole M > Mf and Mf > Mq > Mew- However, we also want to consider the possibility 
that the black hole decays to lower masses, that is Mf <C Mq- We therefore use the phase space 
measure in the Planckian number density for the emitted particles given by 

Nmik) = -^^^^ , (3.1) 

where T denotes the black hole temperature, m the particle mass, k the particle 3-momentum 
and Lo = ^k"^ + m?. Moreover, since Nm depends on m and not just on the statistics, one can 
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no more assume (as in CHARYBDIS) a fixed ratio of production for fermions versus bosons but 
proper particle multiplicity must be used when generating particle types randomly. We used the 
multiplicities predicted by the Standard Model ^32j and modified the code accordingly. 
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Figure 1: Black hole temperature (in GeV) as a function of the black hole mass (in GeV). Left panel: 
comparison between the Hawking law Th (solid line) and Ti from Eq. (j3.3p for I = lO^^GeV^^, 
a = 5, n = 5 (dashed line) and £ = lO'^ GeV-\ a = 5, n = 1 (dotted line). Right panel: Th (solid 
line) compared to T2 from Eq. (fO|) with i = l.li ■ 10~^ GeV~^, ni = 1, 112 = 1 (dotted line). 



In order to include possible quantum gravity effects, we employ modified expressions for the 
temperature T = Ti of the form 



(3.2) 



in which, in order to cover results in the existing literature (for a partial list of approaches to the 
problem, see Refs. [12 \ 114 1 l^lj) we shall consider the functional forms 



Fi{i, a, n) 



T>n 



and 



F2{i,ni,n2) 



1 



Rh 



nil "2 



(3.3) 



(3.4) 



where £, a, n, rii, and 712 are parameters that can be adjusted (see below). Note that -Fi leads 
to a vanishing temperature for vanishing black hole mass (horizon radius) whereas with F2 the 
temperature vanishes at finite M when Rn = i [31j (such a remnant was also considered in Ref. [27j). 
Specific examples which will be used throughout the paper are shown in Fig. [1] together with the 
standard Hawking law (j2.2p . 

A list of some adjustable parameters in the code is given in Table [H The initial black hole 
mass Mq can be either fixed (to within the maximum centre mass energy expected at the LHC) 
or generated according to the partonic cross section. In the simulations we present here, we shall 
take a total centre mass energy of 14 TeV and mostly focus on the the case of Mq = 10 TeV (with 
Mq = 2 TeV taken for comparison 0). At the end point of the decay, the black hole explodes 



^According to the current understanding of black hole production, light black holes of mass around 2 TeV are 
more likely to be produced at the LHC (see Ref. but detailed features in their decay signal would also be less 
clearly identifiable than in those of larger mass. 
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in a selectable number A'f of fragments when M < Aff. We shall use Nf = 2 and set the total 
number of space-time dimensions equal to 6. We shall also set the grey-body factors equal to 1 
for all kinds of particles. This approximation is certainly restrictive and calls for a more refined 
analysis. However, we note that in the standard picture the use of more accurate grey-body factors 
alters the final results only by a small amount and we thus feel that our modified code can be 
regarded as a good tool to study qualitative features related to the late stages of the evaporation. 
A substantial modification of CHARYBDIS will likely be necessary if one wants to obtain more 
accurate predictions about, for example, the actual modification of the Hawking emission. 



Mo 


Initial black hole mass 


1 — 14 TeV; random 


Mf 


Minimum black hole mass 


1 - 1000 GeV 




number of final fragments 


2-6 


d 


number of extra dimensions 


1 - 5 


Fi 


modified temperature 


1 (no modification), Fi, F2 



Table 1: Parameters that can be adjusted in our Monte Carlo generator and their ranges. For each 
modified temperature, one can also set the corresponding parameters as described in the text. 



4 Simulation Results 

The standard CHARYBDIS generator simulates the evaporation of a micro-black hole according 
to Hawking's law until the black hole mass M reaches the minimum value Mf which is assumed 
equal to the fundamental scale of gravity Mq in a given model of extra dimensions {e.g., at 1 TeV). 
The black hole subsequently decays into a few bodies (Af = 1, ... ,6) simply according to phase 
space. As described in the previous section, in our modified code, we implemented ways to extend 
the evaporation to a mass Mf well below the fundamental scale {e.g., down to 1 GeV). 

In the following, we shall compare the outputs of the standard and modified generators from a 
phenomenological point of view. In particular, we shall consider the output of the standard code 
with Mf = 1 TeV (and Mf = 1 GeV for completeness) as reference. The results for the modified 
temperature Tf were obtained with Mf = 1 GeV and the parameters £ = 10~^ GeV~^, a = 5, n = 5 
(see Eq. (13. 3p ) and those for the temperature T2 with £ = 1.14 • 10"'* GeV~^, ni = 1, n2 = 1 (see 
Eq. (j3.4p ). Note that the latter choice of parameters also corresponds to a final mass of the black 
hole equal to Mf ~ 1 GeV at which the temperature T2 = 0. These cases are the same that were 
shown previously in Fig. [TJ 

4.1 Primary Emission 

Let us first examine the "primary" black hole emission, namely the particles produced by direct 
black hole evaporation before parton evolution and hadronization are taken into account. As we 
mentioned before, we set the mass for the final decay of the black hole at Mf = 1 TeV in the 
standard case and at Mf = 1 GeV for the modified temperatures. 

A comparison of the relative abundance of the Standard Model particles produced by the black 
hole evaporation with the Hawking temperature and the modified temperature Ti from Eq. (|3.3p 
is shown in the left panel of Fig. [2j For the temperature Ti, this plot clearly shows a much larger 
multiplicity of isotropically emitted light particles (i.e., photons and neutrinos as well as electrons 
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Figure 2: Left panel: abundance of Standard Model particles produced during the evaporation of 
10 TeV black holes before hadronization, for a run of 10000 events with the temperature Ti (solid 
histogram) and with the standard CHARYBDIS generator (dashed histogram). The fundamental 
scale of gravity is set to 1 TeV, temperature parameters are the same as in Fig. [T] and number of 
extra dimensions d = 2. Right panel: energy distribution of the emitted particles obtained from Ti 
(solid histogram) , including only particles emitted when the black hole has a mass exceeding 1 TeV 
(dashed histogram), and with the standard CHARYBDIS generator (dotted histogram). 

and muons). In the right panel of Fig. [21 we also show the energy distribution of all the emitted 
particles. The modified law dramatically changes the spectrum at low energy (there are about 
three orders of magnitude more particles with energy below 100 GeV than in the standard case), 
leaving the spectrum at large energy moderately affected (by just about one order of magnitude). 
Note that the low energy distribution of particles emitted when the black hole mass M > 1 TeV 
follows closely the standard curve but that there are already missing particles in the high energy 
tail (above ITeV). 

For completeness, in Fig. [3l we show the particle abundance obtained with the standard 
CHARYBDIS generator if we set Mf = 1 GeV and compare it with the results from the modi- 
fied code for the two different temperatures Ti and T2 and Mf = 1 GeV. It then appears that the 
temperature Ti leads to a number of light particles significantly larger than the other three cases. 
The temperature T2 instead produces about the same amount of light particles as the standard 
Hawking temperature provided the evaporation is extended down to Mf = 1 GeV. The correspond- 
ing energy distributions are given in Fig. [H from which one concludes that the number of high 
energy particles produced with T2 is also closer to the standard case with Mf = 1 GeV. To con- 
clude, if one just considers the above data, there is no visible difference in the primary emission 
between T2 and the Hawking temperature provided both are extended to a minimum black hole 
mass of 1 GeV, whereas the spectrum generated with Ti still shows appreciable differences. 

These behaviors can be understood in the following way: since the temperature Ti tends to zero 
for M = 0, we do expect a continuous emission of increasingly softer particles until the black hole 
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evaporates completely (i.e., it reaches the minimum value Mf = 1 GeV). Moreover, once the black 
hole mass has decreased below the threshold for producing a given massive particle, such a particle 
can no longer be emitted. Hence one also expects that the production of the heavier particles 
(i.e., massive gauge bosons and top quarks) be scarcely affected, while the emission of soft low 
mass (or massless) particles should be largely enhanced. For the temperature the situation is 
similar except that its functional form follows the Hawking law more closely and the distribution of 
particles more closely resembles the standard case provided one allows the black hole to evaporate 
down to the same mass at which T2 vanishes (1 GeV in our samples). 

4.2 Final Output 

The primary output includes quarks and gluons that cannot be detected. We therefore used 
Pythia [33] to model the parton evolution and hadronization of quarks and gluons emitted during 
the evaporation (or not participating to the black hole events). This produced the final aspect of 
the events as would be seen by a detector at the LHC. 

In Fig. [5l we show the final particle abundance for the standard Hawking temperature and 
Mf = 1 TeV or Mf = 1 GeV and the modified temperatures Ti and T2 with Mf = 1 GeV. The first 
thing one notes is that the standard Hawking temperature for both Mf = 1 TeV or Mf = 1 GeV 
and the temperature T2 now produce very similar numbers of particles. In contrast to the primary 
output, the temperature Ti now produces a number of photons and pions smaller than the other 
cases. The number of muons and electrons in the final state is however still higher (by a factor 
around 2) for Ti with respect to all other cases. In the case of a black hole with initial mass of 
2 TeV, Fig. [6] shows no evident difference in the number of produced particles for the four cases 
previously discussed. 

Another remarking feature of the events produced with the temperature Ti is that the distribu- 
tion of the sum of transverse momenta for particles with transverse momentum higher than 30 GeV 
is peaked around a larger value (about 5 TeV) compared to the results for the other three cases, 
as shown in the left graph of Fig. [7] (in the right graph, the analogue distributions for Mq = 2 TeV 
are showed). If one counts all particles, however, such distributions show no appreciable differ- 
ence. Finally, we remark that the distribution for the sum of the muon transverse momenta for the 
temperature Ti does not show a peak around zero, contrary to all other cases. The same feature 
appears if one only counts muons with transverse momentum over 30 GeV see Fig. El 

5 Conclusions 

We have considered possible modifications to the standard picture of the decay of micro-black holes 
that might be produce at the LHC. Inspired by the microcanonical description of the Hawking 
evaporation and other theoretical approaches in the literature, we have studied modified statistical 
laws for the emitted particles described by different mass-dependent temperatures Ti and T2 which 
have a regular behavior for vanishing black hole mass. We have also required that energy be 
conserved all the way to the end of the evaporation at a scale Mf <C Mq. We have not explicitly 
considered the case in which black holes leave stable remnants, since that has already been reported 
in Ref. [27]. Further, our temperature T2 corresponds to a zero-temperature remnant of finite mass 
which suddenly decays in a few particles. 

The numerical simulations we have run so far suggest that this kind of modifications, if present, 
might have significant and testable effects for black holes produced with a (relatively) large initial 
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mass (of the order of ten times the fundamental scale of gravity Mq). On assuming that Mq = 
1 TeV, the LHC at full luminosity is expected to produce around three 10 TeV black holes per 
day d]. Such a production rate, although smaller than the 1 Hz for 5 TeV black holes (and even 
larger for 2 TeV black holes) is however large enough to collect sufficient data to compare with 
the results we have presented here and their signature should also be easier to identify than that 
for lighter black holes. For example, one expects to see differences in the number of particles 
produced and in their total transverse momentum. It is therefore important to take into account 
these possibilities in any detailed simulation of the evaporation process whose outputs will have to 
be confronted with forthcoming data. 
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Figure 3: Abundance of Standard Model particles produced during the evaporation of lOTeV 
black holes for a run of 10000 events before hadronization. Upper graphs: standard CHARYBDIS 
generator with Mf = ITeV (left panel) and Mf = 1 GeV (right panel). Lower graphs: modified 
code with the temperature Ti (left panel) and T2 (right panel) and Mf = 1 GeV. 
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Figure 4: Energy spectrum of particles produced during the evaporation of 10 TeV black holes for 
a run of 10000 events before hadronization: standard CHARYBDIS generator with Mf = 1 TeV 
(dotted line) and Mf = 1 GeV (dashed line), modified code with the temperature Ti (solid line) 
and T2 (indistiguishable from the dashed). 
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Figure 5: Final abundance of Standard Model particles produced during the evaporation of 10 TeV 
black holes for a run of 10000 events. Upper graphs: standard CHARYBDIS generator with 
Mf = ITeV (left panel) and Mf = 1 GeV (right panel). Lower graphs: modified code with the 
temperature Ti (left panel) and T2 (right panel). 
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Figure 6: Final abundance of Standard Model particles produced during the evaporation of 2 TeV 
black holes for a run of 10000 events. Upper graphs: standard CHARYBDIS generator with 
M{ = 1 GeV Lower graphs: modified code with the temperature Ti (left panel) and T2 (right 
panel) . 
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Figure 7: Statistical distribution of the sum of particle transverse momenta including only particles 
with transverse momentum higher than 30GoV for 10000 events for the decay of lOTeV black 
holes (left graph) and of 2 TeV black holes (right graph) : standard CHARYBDIS generator with 
Mf = 1 TeV (pointed line) and Mf = 1 GeV (dashed line), modified code with the temperature Ti 
(solid line) and T2 (indistinguishable from dashed line). 
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Figure 8: Sum of the muon transverse momenta in the range 0—10 GeV including only particles with 
transverse momentum higher than 30 GeV for the decay of 10 TeV black holes. Left panel: standard 
CHARYBDIS generator with Mf = 1 TeV. Right panel: modified code with the temperature Ti for 
Mf = 1 GeV. 
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